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Owing to the quadratic nature of the theory, Einstein-Gauss-Bonnet gravity gener- 
ically permits two distinct vacuum solutions. One solution (the "Einstein" vacuum) 
has a well defined limit as the Gauss-Bonnet coupling goes to zero, whereas the other 
solution (the "stringy" vacuum) does not. There has been some debate regarding the 
stability of these vacua, most recently from Deser & Tekin who have argued that the 
corresponding black hole solutions have positive mass and as such both vacua are sta- 
ble. Whilst the statement about the mass is correct, we argue that the stringy vacuum 
is still perturbatively unstable. Simply put, the stringy vacuum suffers from a ghost- 
like instability that is not excited by the spherically symmetric black hole, but would 
be excited by any source likely to emit gravitational waves, such as a binary system. 
This result is reliable except in the strongly coupled regime close to the Chern-Simons 
limit, when the two vacua are almost degenerate. In this regime, we study instan- 
ton transitions between branches via bubble nucleation, and calculate the nucleation 
probability. This demonstrates that there is large mixing between the vacua, so that 
neither of them can accurately describe the true quantum vacuum. We also present 
a new gravitational instanton describing black hole pair production in de Sitter space 
on the Einstein branch, which is preferred to the usual Nariai instantons and is not 
present in pure Einstein gravity. 
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1 Introduction 

In four dimensions, the Einstein-Hilbert action, supplemented with an arbitrary cosmological 
constant, is the unique gravitational action giving rise to field equations that only depend 
on the metric and its first two derivatives pQ. This is no longer true in D > 4 dimensions, 
when the Einstein-Hilbert action is merely a special case of the Lovelock action, satisfying 
the same property [2]. In D dimensions, the Lovelock action is given by a linear combination 
of [^r 1 ] dimensionally extended Euler characteristics, 

[Vl f 

^Lovelock = 2, a nSn, S n = 2~" / d X^f^g S^'''^ .R 0,1 X Cldl . . . FT™ n c „ dn (1) 

n=0 J M 



For any n, S n is a D-iorm that corresponds to a dimensionally extended Euler characteristic, 
it is a topological invariant in D = 2n (trivial for D < 2n), but becomes dynamical for 
D > 2n, and in just the right way as to keep the field equations down to second order. 
This is because the higher order derivatives that inevitably appear under metric variation 
of (CQ) can all be written as total divergences. Lovelock gravity is dynamically different from 
Einstein gravity for D > 4, even though the theories are indistinguishable in four dimensions 
or less (for a review see [3]). 

Clearly Lovelock theory is the non-trivial and unique classical generalisation of GR in 
higher dimensions and this makes it interesting to study in its own right. Not surprisingly it 
has also some relevance to string theory. Higher order curvature terms are in fact known to 
appear in the slope expansion of the heterotic string, generically introducing higher deriva- 
tives in the metric |4J. This will inevitably lead to perturbative ghosts, which are known to 
be absent in the full non perturbative version of string theory. One can avoid this appar- 
ent inconsistency (to leading order) if the slope expansion takes the form of the Lovelock 
action [5]. Indeed, the second order curvature terms are known to take the form of the 
Gauss-Bonnet combination [B|. Because this corresponds to the non-trivial n = 2 Lovelock 
term in ([IT) , we are guaranteed that there are no higher derivatives in the effective string 
action, and therefore no ghosts to second order. 

This has motivated a great deal of interest in Lovelock gravity, and Einstein-Gauss- 
Bonnet (EGB) gravity in particular (see, for example, [Zl El El OH El EEl EH]). This is the 
simplest of Lovelock's extensions of General Relativity, and is described by the following 
action 

]D _, / ~i o A i D i „. r \ T _ d2 a d D«i i D r>abcd /o\ 



Sgb = k d u Xv ^(-2A + R + aL GB ), L GB = R 2 - AR ab R ab + R abcd R a 
Jm 

where k = l/167rG. The action ([2]) is the sum of the first three Lovelock forms, the first two 
giving rise to the well known Einstein-Hilbert term and cosmological constant, the third to 
the Gauss-Bonnet invariant, ^ '—gLcB- The latter term is topological in four dimensions, and 
becomes dynamical in D > 4 dimensions. EGB gravity has been extensively studied in the 
braneworld context owing to a shared "stringy" and cosmological motivation [T4"l [T5l [j~6l [T7] 
and as a modified "scalar-tensor" theory of gravity in 4 dimensions (see for example [T8]). 



In this paper, we will restrict attention to D dimensional EGB gravity, coupled to ordi- 
nary matter fields. The field equations are given by 

G ab + Ag ab + aH ab = —T ab (3) 

Ik 

where T ab is the matter energy-momentum tensor, G a b is the Einstein tensor for the metric 
9 ah, and 

Hgh — 2-Ri? a fe — AR ac R c b — AR acbd R c + 2R acde R b c e — -L GB g ab (4) 

is the second order Lovelock tensor [2]. Given the higher order nature of the theory ([I]), 
the vacuum field equations generically admit multiple solutions with maximally symmetry 
(up to the order n in number). For EGB, we therfore have up to two distinct maximally 
symmetric vacuum solutions, with two possible effective cosmological constants, 



A£ = Acs [l ± ^1 - ^j , (5) 

where 

1 (D-l)(D-2) 

AcS ~ ~4^(D-3)(D-4) (6) 

This is true even in the absence of a bare cosmological contant in the action. For these vacua 
to be well defined, the bare cosmological constant must satisfy the bound, A/A C s < 1/2. It 
is easy to check that A^ s /Acs > 1 > A~ s /Acs, with equality when the bound is saturated, 
A = Acs/2- This is known as the Chern-Simons limit (at least in odd dimensions) [34| . and 
corresponds to the case where the two roots actually coincide. Note that only the lower root 
has a smooth limit, A~ s — > A, as a — > 0, and as such is often refered to as the "Einstein" 
branch. In contrast, the upper root, A^ has no smooth limit as a — > 0, and represents a 
distinct new feature of EGB gravity that is completely absent in higher dimensional General 
Relativity. For this reason, this branch is often refered to as the "stringy" or "Gauss-Bonnet" 
branch. 

Our main interest here lies in the physical relevance of these vacua and hence in their 
stability, in particular that of the Gauss-Bonnet vacuum. There has been much to-ing and 
fro-ing in the literature over the years regarding this issue pj EEl E] • We will review some 
of the various arguments in the next section, before clarifying the situation, we hope, once 
and for all. As a brief taster, note that whilst we agree with Deser and Tekin's recent claims 
that the gravitational mass of spherically symmetric solutions always takes the same sign 
as the inertial mass of the source [UJ , we strongly disagree with their conclusions regarding 
the classical stability of the vacua. We will explain in detail why the Gauss-Bonnet vacuum 
is certainly itrTstable, at least perturbatively, and why this does not conflict with Deser 
and Tekin's calculations of gravitational mass. In short, the instability is due to a freely 
propagating tensor ghost that is not excited by the spherically symmetric solutions. 

It turns out that the perturbative analysis of section [2] actually breaks down close to the 
Chern-Simons limit, due to strong coupling. At this point, we need to investigate Gauss- 
Bonnet vacuum stability using non-perturbative techniques. To this end, in section [31 we 



will focus on instanton solutions, corresponding to differing vacuum bubbles surrounded by 
domain walls of ordinary matter. Happily, our description works best in the near Chern- 
Simons limit, in contrast to the analysis of section [2J By employing standard semi-classical 
techniques [201 [21], we will find that the mixing between branches is completely unsuppressed. 
Our results indicate that in the near Chern-Simons limit, an empty, homogeneous vacuum, 
on either branch, does not describe a stable quantum vacuum, as it becomes littered with 
bubbles of vacua from the other branch. 

We conclude this paper with some discussion, presenting a novel everywhere smooth 
EGB gravitational instanton with a non-zero temperature. This solution exists only on the 
Einstein branch, but is absent in General Relativity. It can be used to demonstrate the 
semi-classical instability of de Sitter space, in addition to the standard Nariai instanton. 

2 Chronicles of a stringy vacuum 

Let us begin this section with a short account regarding the status of Gauss-Bonnet vacua 
and their physical relevance. We begin, not with the vacuum itself, but with a spherically 
symmetric source of inertial mass, M. Now in Einstein gravity, we know that Birkhoff's 
theorem applies outside the source, and that the corresponding solution is Schwarzschild or 
Schwarzschild-(anti) de Sitter 

dv 
ds 2 = -V{r)dt 2 + — - + r 2 dQ D _ 2 (7) 



V(r) 



where the potential is given by 



2Ar 2 1 M 

V(r) = — — — — + 1 



(D-l)(D-2) (D - 2)Q D _ 2 Kr D - 3 v/ 

and dVt£>_2 is the metric on a unit (D — 2)-sphere. The gravitational (ADM) mass, E of this 
solution exactly matches the inertial mass of the source, i.e. E = M [22j. As long as M > 0, 
the singularity at r = is shielded by a horizon in accordance with cosmic censorship q In 
Gauss-Bonnet gravity, a generalised form of Birkhoff's theoreno also holds [HJ HSl [23], and, 
as expected, there are two branches of Schwarzschild-like solutions, with potentials given 

by mi 



T , / x 2A C5 r 2 / / 2A M(D-l) \ 

mo = 1 - P _ x)( g _ 2) (' * f - a^ - A JJ- J ■ < 9 > 

Now, even for M > 0, we immediately see that the cosmic censorship hypothesis will not 
always hold on the Gauss-Bonnet branch - naked singularities can form even for sources 



^n D = 5 there is critical bound for the mass beyond which the solution is singular but this does not 
affect the argument presented here 

2 In the Chern-Simons limit, an extra class of degenerate solutions exist [TH] 



of positive inertial mass! The situation becomes even more bizarre when one considers the 
asymptotic behaviour of the potential §B) 

2A^ ff r 2 M , s 

V+(r) ~ — h 1 ± h (10) 

One is tempted to compare this with the potential (jSJ) in Einstein gravity, and to naively 
conclude that the gravitational mass is given by 

E± = ^M (11) 

If this were correct, it would mean the solution on the Gauss-Bonnet branch had negative 
gravitational mass, even though the source had positive inertial mass! Again, if correct, this 
would certainly render the Gauss-Bonnet vacuum unstable, as it would not be a true local 
ground state of the theory. This weird behaviour was attributed to perturbations around 
the vacuum containing a ghost [7J [H [10]. To see the ghost, consider metric perturbations, 
$9ab = h ab about the maximally symmetric vacuum solution, g ab , with effective cosmo logical 
constant, A eff . Expanding the action (|2]) to quadratic order, we find that 



6 2 S = --k[1--^) / d D xV^h ab (S 1 G ab + A efi h ab )+6 2 S m (12) 




where 5iG ab is the Einstein tensor expanded to linear order in h ab , and (^m is the quadratic 
contribution from the matter part of the action. Variation of this action gives the linearised 
field equations 

1 - ^-) (S.Ga, + A cS h ab ) = ^5 x T ab (13) 

Acs J 2k 

The Einstein branch is well behaved at this order, since A e s/Acs < 1> and so the kinetic part 
of the gravitational action has the usual sign. Indeed, we see from the equations of motion 
( Tl3l) that the linearised theory is the same as in Einstein gravity, with a positive effective 
Newton's constant. In contrast, on the Gauss-Bonnet branch, A e g/Acs > 1> so the kinetic 
term in the action has the opposite sign, indicating the presence of a ghost. We also see 
how the linearised theory corresponds to Einstein gravity with a negative effective Newton's 
constant. This corresponds to an antz-gravity theory and leads to gravitational repulsion, 
consistent with the expansion of the potential fllOp . 

Because the kinetic term of a ghost has the wrong sign, we can think of it as carrying 
negative kinetic energy. As a result, it was originally claimed [7J [HI [10] that the ghost was 
responsible for the negative gravitational energy stored in the gravitational field around a 
spherically symmetric source of positive mass, as indicated by equation ( ITTj) . However, this 
statement cannot be correct, at least perturbatively around the vacuum, as we will now 
explain. 

A ghost can only lower the gravitational energy of the gravitational field if it is freely 
propagating. Now, the perturbation about the Gauss-Bonnet vacuum can be split into 
tensors, vectors and scalars in the usual way. By direct comparison with Einstein gravity, 



which is equivalent up to a sign at linear order, we can conclude that only the tensor modes 
propagate freely on the vacuum. Put another way: there are no scalar or vector gravitational 
waves. So, for the energy of the gravitational field to be negative, the spherically symmetric 
solution (Q must excite the negative energy gravitational waves coming from the vacuum 
tensor perturbations. However, this is not the case. The spherically symmetric solution only 
excites scalar modes. Since these do not propagate freely, by virtue of Birkhoff's theorem, 
they only carry the energy given to them by the source. The gravitational energy of the 
black hole can only depend on the overall sign of the gravitational Hamiltonian, and not on 
which branch you are on. 

This statement is actually consistent with some of Deser and Tekin's more recent work |llj . 
By developing new techniques for calculating the ADM and AD energy of quadratic gravity 
theories, they found that the correct energy expressions should actually be given by 

E± = M, (14) 

rather than equation ( fTTT) . The gravitational energy is always positive for a source of positive 
inertial mass, whatever branch you happen to be on! These new expressions were later 
confirmed using the Gauss-Bonnet Hamiltonian [12], and agree with the reasoning given in 
the previous paragraph. 

Deser and Tekin went on however, to claim that the positivity of the energy (IT4"1) meant 
that the vacua on both branches were stable. We do not agree with this claim. We have 
already explained how the spherically symmetric solutions only excite scalars. These do not 
propagate freely and can only carry the energy given to them by the source itself. However, 
we cannot rule out the possibility that other, non-spherically symmetric, solutions exist that 
excite the freely propagating tensor modes, eg a binary system with a time dependent mass 
quadrupole moment. On the Gauss-Bonnet branch, these do indeed give rise to negative 
energy gravitational waves, rendering the Gauss-Bonnet background unstable. 

This kind of instability can manifest itself already at a classical level. To see how this 
might occur consider, for example, a binary system of black holes, or indeed any other source 
capable of emitting gravitational waves in D dimensions (see [24J for some recent discussion 
of gravitational waves and binary systems in higher dimensional GR). If our theory of gravity 
happens to contain freely propagating tensor ghosts, the energy of our source will increase 
as it emits waves carrying negative energy. For the case of the binary system, its orbital 
period will become quicker and quicker as more and more gravitational waves are emitted 
at an ever increasing rate! 

Even in a true vacuum, with no matter excitations, a ghost will lead to a quantum 
instability if it is coupled to matter, as it will be spontaneously produced in the vacuum 
along with ordinary matter fields [22]. Since the ghost carries negative energy, and the matter 
fields, positive energy, this can occur without any overall violation of energy conservation. 
For a Lorentz invariant theory with no UV cut-off, the rate of particle production diverges, 
and the vacuum is destroyed infinitely quickly (see [26]). 

Our conclusion, then, is that although the Einstein vacuum is stable, the Gauss-Bonnet 
vacuum is certainly unstable. This perturbative instability is due to a tensor ghost, which 
can propagate freely on the vacuum. It can be seen even at a classical level whenever there 
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is a source of gravitational waves. In the absence of any excited sources, the vacuum is still 
quantum mechanically unstable as long as the ghost couples to ordinary matter fields. Note 
that the ghost is not excited by the spherically symmetric solutions given by (J7J and ([HD, 
so it comes as no surprise that Deser and Tekin ultimately found the gravitational energy 
of these solutions to be positive. Those solutions only excite scalars, which are not freely 
propagating, and therefore only carry energy given to them by the source. 

In the interests of completeness, we note that all of these results rely on the fact that the 
overall gravitational coupling, k = 1/16tcG > 0. If we reverse the sign of k (without altering 
the matter lagrangian), the situation is reversed: the Einstein branch suffers from a ghost, 
whilst the GB branch is stable. Since the overall sign of the gravitational Hamiltonian is 
also reversed, the black holes masses take the opposite sign to the inertial mass of the source, 
but again, do not depend on which branch you are on. This theory does not admit a sensible 
general relativistic limit as a — > 0, as the Einstein-Hilbert piece always takes the wrong sign 
in the action, so we will not discuss it further. 

We end this section with a note of caution. The perturbative analysis we have presented 
here is not universally applicable. In particular, we see that it will fall victim to strong 
coupling when the coefficient of the kinetic term in ( Tl2l) becomes very small, i. e. as A e ff — > 
Acs- This corresponds to nearing the Chern-Simons limit, when higher order interaction 
terms become important, and one needs to go beyond perturbation theory in order to study 
vacuum stability. This will be the subject of the next section. 

3 Bubbling Gauss-Bonnet vacua 

In the last section, we saw how the Gauss-Bonnet vacuum suffered from a perturbative ghost 
instability. This result holds firm as long as we are not too close to the Chern-Simons limit, 
at which point perturbation theory breaks down due to strong coupling. In this section, 
we will study vacuum stability using non-perturbative phenomena. We will be interested 
in instantons that represent tunnelling between vacua, through bubble nucleation. As we 
will see later, our description will actually work best close to the Chern-Simons limit, i.e. 
precisely when the analysis of the previous section becomes invalid. The techniques for 
studying vacuum decay were developed by Coleman and his collaborators, first for QFT in 
Minkowski space [20] , and later including the effect of four dimensional Einstein gravity [21] . 
We will now extend that analysis to D dimensional Gauss-Bonnet gravity (see also the very 
recent study [27]). 

We begin by constructing the so-called "bounce" geometries, named after the correspond- 
ing entity in particle mechanics. Typically, these describe spherically symmetric bubbles of 
"true" vacuum expanding in the "false" vacuum. Occasionally, the reverse is possible, with 
bubbles of false vacuum forming inside the true vacuum with a small but non-zero probabil- 
ity [SO]. In the thin wall limit, which we will consider here, the two vacua are separated by 
a domain wall, composed of ordinary matter. 

Let us denote the interior of the bubble by M.\ and the exterior by M.2- We will assume 



that for i — 1, 2, M.i is a maximally symmetric solution to the vacuum field equations 

G ab + Aig ab + aiH ab = 0. (15) 

Note that at this stage we are not assuming that Ai = A 2 or a,\ = a.%. We do this for 
the purpose of generality, since we might want think of Gauss-Bonnet gravity as an effective 
theory, with these quantities corresponding to the vevs of other fields. However, our ultimate 
interest will lie in transitions between branches, and the stability of the vacua for a given 
set of vevs. So, although we will keep things general for the most part, we will occasionally 
assume that the bare parameters are indeed the same on either side of the wall. 

The wall itself represents the common boundary of the two manifolds, and is denoted 
by S = dM.\ = dM.2- We will assume that the wall has tension a, so that the junction 
conditions^ there are given by [28J 



kA 



K^v - K^ pv + 2a ( Q^ - -Q-y^ 



h^ ( 16 ) 



where the jump operator is AX = X 2 — X\ and the averaging operator is (X) = (X 1 + X 2 )/2. 
The extrinsic curvature of the wall in Jv[ { is given by the Lie derivative of the induced metric 
7^, with respect to the unit normal, n a , to the wall pointing out of M.i and into M. 2l 

K$ = \tnl» V . (17) 

If the Riemann tensor on the wall is given by H^ap, then, suppressing the label (i), 

Q^ = 2KK flX K\ + K Xp K xp K llu -2K llX K Xp K l/p -K 2 K^ 

+2KK^ + KK^ - 2K xp TZ Xlipu - ATZ^K^ (18) 

Note that we have used Latin indices to describe the bulk coordinates, and greek indices to 
describe the wall coordinates. In particular, it is convenient to introduce Gaussian-Normal 
coordinates x a = (£,x M ) relative to the wall, which is fixed at £ = £o- The bulk metric is 
then given by 

ds 2 = g ab dx a dx b = d£ 2 + p(£) 2 [-dr 2 + cosh 2 rdVt D _ 2 ] (19) 

Note that the expanding bubble wall corresponds to a de Sitter hyperboloid of constant 
radius, p(£o)> embedded in the bulk geometry. If the effective cosmological constant in Aii 
is written as A^ = —\{D — l)(D — 2)kf, then the function p(£) is given by 

p(0 = 7 ^sinh(fc 1 0, 0<£<£ (20) 



3 Junction conditions are uniquely determined in Lovelock theory due to the fact that equations remain 
second order and boundary terms to a manifold with boundary are uniquely defined as dimensionally ex- 
tended Chern D — 1 forms (see for example 3.). The relevant boundary term for Gauss-Bonnet was given 
by Myers in [Sj. Varying this boundary term gives the result of [5S] 



in A^i and 



in 7M2, where 



p(0 = r sinh[fc 2 (£ -/?)], &<£<& 

&2 



& 



OG 



for k\ > 



/5 + vr/|A; 2 | forA;2<0 



f21^ 



f22) 



Note that these formulas are valid, even when kf < 0, as long as the appropriate limit 
or analytic continuation is understood. The boundary conditions at the wall are given by 
continuity of the induced metric 



A 



sinh kX 

k 



0, Ax = 6, A 2 = £o-/3 



(23) 



and the junction conditions (TT6I) . The latter gives the following expression for the wall 
tension 



a = -2(D-2)kA 



P -\l-\{D-Z){D-A)a 



P_ 

P J P z 



(24) 



e=?o 



This expression only really makes sense if the thin wall approximation is a good one. For 
this to be so, we require [2TJ 



A//(£o) 



p(&) 



< 



1 



Ptfo 



A[coshfcA] 



<1 



&2U2 - IfcilAi 



<1 



(25) 



Furthermore, the wall is assumed to be made up of ordinary fields, and therefore its energy 
momentum tensor, T ab = —8(£ — ^crj^S^, must satisfy the weak energy condition. This 
means that we require o > for viable bounce geometries. As a check, consider the case 
where the vacuum is Einstein on both sides of the bubble wall, and take the GR limit, a\ = 
«2 = 0, so that a = —2{D — 2)kA [cosh/cA] /p(£o)- Assuming the thin wall approximation 
(j25p . we see that a > 0, if, and only if, we have one of the following 

(i) Ai < A2 < 

(ii) Ai < < A 2 

(iii) 0< Ai< A 2 , |A;i|Ai< |A; 2 |A 2 < f 

(iv) 0< A 2 <A l5 |A; 2 |A 2 > |fci|Ai > f 



Cases (JTj) to (Jul]) are consistent with [21] . and describe tunnelling between a false vacuum of 
large A, to a true vacuum of smaller A. If both vacua are de Sitter, we see from case (jrv|) 
that tunnelling can occur in the opposite direction, from true to false, but with a very low 
nucleation rate (see, for example [30]). 
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We are primarily interested in transitions between branches in Gauss-Bonnet gravity. 
With this in mind, let us consider the case where the bare parameters are the same inside 
and outside of the bubble (i.e. AA = Aa = 0). It then follows that 

ff = -H(Z)-2)( C -3)(Z)-4) K a(^^) 3 (26) 

The thin wall approximation (1251) has extra significance here, since it corresponds to being 
close to the Chern-Simons limit. This is precisely when the perturbative analysis of the 
previous section was in danger of breaking down due to strong coupling effects. Because 
they have opposite regimes of validity, we see explicitly how the two sections compliment 
one another perfectly. On the one hand we have a perturbative analysis that works best far 
away from the Chern-Simons limit, whereas on the other hand we have a non-perturbative 
analysis that works best close to that limit. 

Given the expression (TSHJ), and assuming the thin wall limit approximation (TSo'j) . we see 
that a > 0, if, and only if, one of the following holds 

(a) a > 0, Af < Af < 

(b) a > 0, Af < < Af 

(c) a< 0, < Af < Af , lifcxlAi > |A; 2 |A 2 > f 

(d) a < 0, Af < < Af 

(e) a < 0, < Af < Af , \k 2 \\ 2 < \ki\Xi < § 

Cases (jaj) to (jcj) describe tunnelling from high to low effective cosmological constant, con- 
sistent with decay of the false vacuum through the nucleation of a bubble of true vacuum. 
The remaining cases, ((dl) and (|ej), describe the reverse process, in which true vacuum decay 
occurs, and we have tunnelling to higher A e g. These are analogous to case (jTvj) in the GR 
limit, although now we no longer require both vacua to be de Sitter. 

Which is the true and false vacuum depends on the sign of a. For a > 0, the Gauss- 
Bonnet vacuum may be thought of as the true vacuum, as we always have A^ < Acs < A^, 
whereas for a < 0, it is the Einstein vacuum that is the true vacuum, as A^~ ff < Acs < ^-ts- 
Note that while case (jcj) corresponds to the nucleation of a bubble of Einstein vacuum in 
Gauss-Bonnet, all other cases correspond to the reverse: bubbles of Gauss-Bonnet vacuum 
in Einstein. 

Vacuum decay occurs when a region of vacuum is replaced with a region of another 
vacuum, such that there is a net energy loss inside the bubble. This energy is used to excite 
the fields that form the bubble wall, so that overall, energy is conserved. We can illustrate this 
explicitly in the present case by calculating the energy of the bubble relative to the region of 
initial vacuum it replaced. This can be done using the Hamiltonian formulae derived in [12] . 
Note that in ffl9l) we have expressed the bulk metric using cosmological coordinates along 
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the direction of the wall. Since we want to calculate an energy, its convenient to transform 
to global de Sitter coordinates along the wall, so that the bulk metric is written as 



ds 2 = d£ 2 + p(£f 



(l-r 2 )dt 2 + 



dr" 1 



1 — -r 



+ r 2 dn D _ 3 



(27) 



Surfaces of constant t, E t , provide a foliation of the bulk with respect to the timelike Killing 
vector d/dt. It follows that the lapse function N = pyl — r 2 , and the shift vector vanishes. 
We also need to define the surface S t , which is the intersection of the wall with E t , and 



has induced metric qijdx l dx^ = p 2 (£o 
inside the bubble is given byj 



dr 2 

1-r 2 



+ r 2 dQ,D-3 ) ■ Now, from [12], the energy stored 



E 



bubble 



d D ~ 2 x ACt 



St 



bdy 



where, for vanishing shift vector, 



C 



bdy 



VqN 



2/C + 12a6 [ ;S?6? K\ (W k mn {q) - 27&W*) - -K\KtK* 



f28) 



(29) 



Here lZijki(q) is the Riemann curvature on St, and /Qj and Tiij are the extrinsic curvatures 
of St in E 4 and the wall respectively. After a little algebra and integration, one finds that 
the energy inside the bubble is given by 



E 



bubble 



vD-1 



A 



2nQ D _ 3 p(£ 

rr 



l--(D-3)(D-A)a 



C=Co 



Q, -"= ^- J 



(30) 



(31) 



where we have used equation ff24l) . When the tension is positive, the energy of the bubble is 
negative relative to the region of vacuum it replaced. This applies in each of the five cases 
(m) to (|ej), regardless of whether or not the bubble is one of true or false vacuum. It couldn't 
be any other way, since the energy stored in the wall must be compensated for by an energy 
deficit inside the bubble. In fact, 



E, 



wall 



St 



d u - l x^dNa = -E } 



bubble 



(32) 



and so the net energy is indeed zero. 

We will now apply standard semi-classical techniques [201 IZE| in order to derive the 
probability of bubble nucleation in the general case. This is given by 



V oc e- B ' h 



(33) 



4 The alert reader will notice that there is sign difference between equation (|28|) . and the corresponding 
energy formula (98) in [12]. In [12] . for any quantity X in the test spacetime, with corresponding value X 



in the background, AX = X — X. In this paper, however, the bubble is given by .Mi, and so X 
X = X2, which means our definition of A differs by a sign, compared with [12] . 



Xi and 
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where the instanton action, B, is given by the difference between the Euclidean bounce 
action and background action 



B = S } 



bounce 



S ] 



background 



(34) 



A general bounce is made up of the bubble interior, .A/f i, the exterior, A4 2 , and the bubble 
wall, E. In Lorentzian signature it is described by equations ( fi~9l to ([211 . To go to Euclidean 
signature we Wick rotate the time coordinate, r — > ir E , so that the Euclidean bounce 
geometry is given by 

ds 2 = g ab dx a dx b = d£ 2 + p(£) 2 ^D-i (35) 

As usual, the wall, which was a de Sitter hyperboloid in Lorentzian signature, corresponds 
to a (D — l)-sphere in Euclidean signature. Taking care to include the relevant surface terms 
in the action [291 E] > we find that the Euclidean bounce action is given by 



S ] 



bounce 



-« / V9 (R - 2A + «£gb) + / y/l k [2AK + \ A{aQ) ] + 

Jmi\jm2 Jz L V 6 J 

+ (boundary terms at infinity) 



a 



(36) 



To describe the background geometry, Ai, we simply extend the solution in Ai 2 all the way 
to the centre of the space, so that the Euclidean background metric is given by 



ds 2 = g ab dx a dx b = d£ 2 + p(^) 2 dfi D _i 



where 



p(0 



k, 



sinh[/c 2 (^ — 



P < z < & 



(37) 
(38) 



The background action is therefore given by 



S ] 



background 



~ K / Vd (-R — 2A + «Lgb) + (boundary terms at infinity) (39) 

JM 



It is clear that when we calculate the instanton action, B, the boundary terms at infinity 
cancel one another. Plugging everything in, and making use of the expression ( 1241) . we find 
that the instanton action is given by 



B 



k£Id-iA 



(-2A - D(D - l)k 2 + D...(D- 3)ak 4 ) t drj ( S J^^ 



+2p u - 2 p <l-2(D-2)(D- 3)a 



P l 



(40) 



C=Co 



This expression is valid even when AA 7^ 0, and/or Aa 7^ 0. Again, we use the GR 
limit as a check, taking both vacua to be Einstein and setting ol\ = a 2 = so that B = 
2k[D — 2)f^D_iA J Q drj{k~ 1 sinh ki]) D ~ 3 . In the thin wall approximation (1251) . we set Aj = 
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(A)(l + 6j), where |ej| <C 1, and then it is straightforward to check that (e) = 0. From (!20l) 
and (12T|) . we see that ki\ = sinh" 1 kip(£ ), where p(£o) is fixed, and so 

B = -^jKn D ^p(^) D /«fcA»AA + 0(e 2 ) (41) 

where f{z) = -— ^ (-^ - ^T Jo dw (f^r)^ 3 )- Now for x, yeRwe have /(a;) > 
for x > 0, whereas /(ii/) > for y e (0, 7r/2) and /(zy) < for y e (7r/2, 7r). As expected [2Tj . 
£> > for each case (EJ) to (jivj), so the probability of tunnelling between vacua is always 
exponentially suppressed. True vacuum decay, given by case (By]) , is even further suppressed, 
as can be seen using the identity f(iy) + f(i{i: — y)) = — sinD J^ dw sm D ~ 3 w [30]. 

Returning to Gauss-Bonnet gravity, we are mainly interested in transitions between 
branches, and the tunnelling processes described by (jgj) to (jej). Explicity setting A A = 
Aa = 0, we find that the instanton action simplifies, a little, to give 




D-l 



(-2A - D(D - l)k 2 + D...(D- 3)ak 4 ) t d V f 8 J^!l\ 



D-l\A[coskkX] t{D _ 2){D _ 3) JM«*W :i 



D-4J pfo) y /v ' \ p(6 



^42) 



As we have already stated, the thin wall approximation can only really be trusted close 
to the Chern-Simons limit, i.e. when A^ = Acs(l + €»), where |ej| ^C 1. For transitions 
between branches the average effective cosmological constant, (A e ff) = Acs, and so (e) = 0. 
For cases (jgj) to (jgj), it follows from ([H that a = C(e 3 ) > 0, and from flU) that B = C(e 3 ). 

Working to order e, we now see how a tensionless domain wall mediates transitions 
between branches in Gauss-Bonnet gravity, altering the effective cosmological constant. At 
this order, the transition probability (|33|) is not suppressed, in contrast to similar transitions 
in General Relativity. This is the main result of this section, and can be interpreted as 
follows: close to the Chern-Simons limit, there is large mixing between the two distinct 
vacua in Gauss-Bonnet gravity. Neither the empty Einstein vacuum, nor the empty Gauss- 
Bonnet vacuum provide a good description of the stable quantum vacuum, as they will both 
become littered with bubbles of the other. 

It is interesting to note that a similar behaviour was seen in the DGP model [31] . In that 
case there is large mixing with empty five-dimensional space and one filled with so-called 
self-accelerating branes [32]. This pathology stemmed from the presence of perturbative 
ghosts on the self-accelerating branch [33J , so given the results of the previous section, it is 
tempting to think that something similar is happening here. This may be so, but there are 
a number of reasons to tread carefully before making such a conclusion. Firstly, as we saw 
in the last section, the perturbative ghost only haunts the Gauss-Bonnet branch, whereas 
here, tunnelling from both branches is unsuppressed. Secondly, and most importantly, this 
section and the last have opposite regimes of validity, and in particular, one cannot really 
trust the existence of the ghost beyond the strong coupling scale, which is very low close to 
the Chern-Simons limit. 
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4 Discussion 

Lovelock theory in D > 4 spacetime dimensions differs dynamically from the naive extension 
of General Relativity in higher dimensions. Unlike GR and Chern-Simons gravity [31] , where 
the vacua are unique, Lovelock theory has multiple vacua for a given set of bare parameters 
in the action. In the simplest non-trivial extension of Lovelock theory, namely Einstein- 
Gauss-Bonnet (EGB) gravity, there are two such vacua: one with a well defined limit as 
a — > (the Einstein vacuum), and one without (the Gauss-Bonnet or "stringy" vacuum ). 
We have demonstrated that the Gauss-Bonnet vacuum is perturbatively unstable, despite 
the fact that the corresponding black hole excitations have positive mass, as shown recently 
by Deser and Tekin [TT]. The instability is due to a perturbative spin 2 ghost that is 
not excited in spherically symmetric vacua owing to an analogue of Birkhoff's theorem in 
Lovelock gravity [T4"l [T5l 123] . However, since the ghost couples to matter, and other ordinary 
fields, it will inevitably lead to a quantum instability of the vacuum through ghost-nonghost 
pair production, without violating conservation of energy. 

At the classical level, the pure vacuum itself is not unstable. A classical instability can 
appear if a small amount of matter is put in by hand so that you are not quite in vacuum. 
To see this, consider a matter source that is capable of emitting gravitational waves such as 
a binary system of black holes with a time dependent mass quadrupole. Since the gravitons 
are ghost-like they carry negative kinetic energy, and their emission leads to an increase in 
the energy of the source. We soon end up with a spacetime containing a very high energy 
source surrounded by a turbulent sea of gravitational waves, and as such, this can no longer 
be regarded as a perturbation about the vacuum. In this sense the introduction of a small 
matter source has rapidly destroyed the vacuum, destabilising the system already at the 
classical level. 

It is natural to ask what is the end point of the ghost instability? One might specu- 
late that it leads to a transition between branches, which may certainly be possible at the 
quantum level. We saw in section [3] how quantum transitions between branches can occur 
through bubble nucleation, although it was not clear whether or not this had anything to 
do with the ghost, or some other superselection mechanism. Another possibility is that the 
ghost mediates vacuum decay into spacetimes of less symmetry, as one might infer from 
the enhanced production of Nariai black holes in six dimensions [17]. We will discuss this 
possibilty in more detail later on. At the classical level, things are just as speculative. Clas- 
sically, it is very difficult to see how a transition between branches could happen since it 
would require a change in the asymptotics. In this case, the end point, if it even exists, 
will presumably depend on the nature of the source, and one may have to resort to complex 
numerical simulations to shed any more light on the issue. 

All these statements regarding the (in)stability of vacua at both the classical and quantum 
level are reliable as long as linearised perturbation theory works well on all but the shortest 
distance scales, i.e. as long as A c s/Acs "C 1. However, as we approach the Chern-Simons 
limit in which the two vacua become degenerate (A e s/Acs ~^ 1), the coupling of the graviton 
to other fields and to itself becomes stronger and stronger. This leads to a breakdown of 
linearised theory at larger and larger scales, pushing out the region of space in which we 
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trust our perturbative analysis. Indeed, in the exact Chern-Simons limit, we see from ( ITBl 
that linearised theory breaks down altogether, on all scales. Close to this limit, we cannot 
make any robust conclusions regarding stability from a perturbative analysis alone. This 
limit has qualitative similarities to the partially massless limit for massive gravity in de 
Sitter [351 [36], aswell as the chiral limit of topologically massive gravitja in AdS [37J 138] - 
and the zero tension limit of DGP self-acceleration [3T| I3"2l [33] . In each case there is a ghost 
whose kinetic term seems to disappear [351 S3 [33] in the appropriate limit and a more careful 
analysis is required to identify any left over degrees of freedom. For massive gravity in de 
Sitter, the dangerous mode genuinely disappears [36], whereas in the latter two cases the 
ghost mixes with another mode, such that one ghost-like degree of freedom remains [38, 33J. 
We would need a direct analysis beyond the scope of this paper to see what happens in the 
exact Chern-Simons limit of EGB gravity. 

Although linearised analysis breaks down close to this limit due to strong coupling, this 
is precisely the regime in which we can reliably study instanton transitions between the 
two distinct vacua arising from the same theory (with the same bare parameters). The 
transitions are achieved through bubbles of one vacuum being nucleated within a sea of the 
other vacuum. The instanton analysis relies on the thin wall approximation, which works 
best near the Chern-Simons limit. This is because the two vacua are energetically very close 
to one another, and so only a small kink is required in the geometry of the bounce. In the thin 
wall/near Chern-Simons limit, we saw that the transition probability of bubble nucleation 
is not (linearly) suppressed unlike similar transitions in Einstein gravity [21] , and concluded 
that there is very strong mixing between the two distinct vacua of Gauss-Bonnet gravity. For 
this reason, in the strong coupling regime, neither vacuum can accurately describe the true 
quantum vacuum state, since both will quickly become littered with bubbles of the other 
vacuum. 

We have focussed on two pathologies associated with vacua in Einstein-Gauss-Bonnet 
gravity: a ghost-like instability in the perturbative regime on the GB branch, and unsup- 
pressed mixing between vacua in the strong coupling regime. A third pathology was recently 
identified on the Einstein branch, through calculations of the viscosity-entropy bound J3U] 
for gauge theories with a Gauss-Bonnet gravity dual. In the background of an EGB planar 
black hole, it was found that graviton wave packets can propagate with superluminal speeds 
whenever ^ < 8ak 2 < 1, where 1/k is the AdS length. This is not necessarily a sign of 
instability of the background, but one might expect some transition and novel physics to 
appear whenever Sak 2 ~ Jr. This scale does not show up in our calculations because we 
have considered perturbations of maximally symmetric backgrounds where the perturbation 
operator f[T3|) is simply a constant multiple of the Einstein perturbation operator. This is 
no longer the case when the background has a non-trivial Weyl tensor, as in the case of the 
planar black hole calculation of [39] . 

We close off the section by complementing our instanton analysis of section [3] with some 
exact gravitational instantons present for positive effective cosmological constant. Some of 
these are present in General Relativity and correspond to the spontaneous creation of black 



We thank Stanley Deser and Bayram Tekin for drawing our attention to this. 



15 



hole pairs in the de Sitter vacuum, which are then accelerated away from one another by the 
background cosmo logical constant. In four dimensional GR, the well known "perdurance" 
calculations [40] demonstrate that de Sitter space is semi-classically unstable to the creation 
of Nariai black hole pairs, albeit at an exponentially suppressed rate (see [31] for the extension 
to higher dimensional GR). In five dimensional EGB gravity, we find that the Nariai black 
hole can only exist on the Einstein branch, and only when the bare cosmo logical constant, 
A = 6H 2 > 0. The corresponding instanton has geometry S 2 x S 3 , and is given by the 
following metric 

^ = i 1 -^) du 2 + -Lm, (43) 



AH 2 J 2H 2 

Recall that A /Acs < 1/2 and so 1 + 8aH 2 > 0. This means that the metric (1331) is well 
defined except in the Chern-Simon limit when the two-sphere shrinks to zero size. The 
background geometry is obtained by Wick rotating the Lorentzian de Sitter vacuum to 
euclidean signature. It is therefore given by a five-sphere of radius H~^ where 



2 A cff VI + 8aH 2 - 1 
HcS = T = Ta (44) 

is the effective curvature of the background. The instanton action ( 1331) is given by 

B Narim = 4nn 3 K I h3 J (45) 

which is always positive in the regime for which the solution (1331) is valid (H 2 > 0, 8aH 2 > 
—1). The creation rate of Nariai black hole pairs is therefore exponentially suppressed, since 
V oc e~ BNariai ^ h , as in General Relativity. 

Interestingly, it turns out that there is an alternative channel for the creation of black 
hole pairs in five dimensional de Sitter, unique to EGB gravity. The process is described by 
a novel gravitational instanton solution that is, again, only present on the Einstein branch, 
and has metric 

ds 2 = V{r)dr 2 + y— + r 2 dQ 2 3 (46) 

with 




V(r) = 1 ! - ( 1 - \l (1 + 8aH 2 )(l + ^) 1 ( J7) 



r i 



The solution is only valid when a > and < H 2 < l/8a, and is smooth everywhere in 
the range, r_ < r < r + where r\ = lj: 2 ~^° ~ are ^ ne horizon positions, V(r±) = 0, for 
the corresponding Lorentzian black hole. The smoothness of the solution (1461) is guaranteed 
by the fact that V'(r + ) = — V'(r_), which means both horizons have the same temperature 
T = |l // ( r ±)l/4 7r - Upon Wick rotating the Lorentzian black hole solution to Euclidean 
signature, this ensures that the instanton solution f|36l) is free from conical singularities at 
both r_ and r + , provided we choose Euclidean time to have period 1/T. The instanton (136"]) 



16 



describes the sponteneous creation of the Lorentzian black holes in the de Sitter vacuum. 
The creation rate is again governed by the instanton action, 

'1 + 12«# e 2 ff 



J^alt.p.rn.nt.iiip. 4r7Ti oSrv 



H 3 



— (r+ + r. 



(4* 



As before, this is always positive within the regime of validity of the solution (1461) . so the 
corresponding black hole creation rate is again exponentially suppressed. However, when 
both instantons solutions, (J4"3"l) and (j4"6| ). exist, it turns out that B^ ariai > B a i ternative , with 
equality^ at the endpoint H 2 = l/8a. This means that the creation of these alternative 
black holes exceeds that of the Nariai black holes! In fact the creation rate of the alternative 
black holes approaches order one as a — > 0. This is because the alternative black hole, which 
has mass M a i ternative = 6a(l + 8aH 2 )Q 3 K, becomes arbitrarily small in the GR limit. 

We wish to reiterate the fact that no such gravitational instantons exist on the Gauss- 
Bonnet branch in D = 5 dimensions. However, in D = 6, as was recently shown in [IT] . 
there exist Gauss-Bonnet branch Nariai instantons for which the instanton action is negative 
and therefore black hole production is exponentially enhanced! Naively, this would seem to 
reflect the presence of the perturbative ghost discussed in section [21 Of course, this ghost is 
present in D = 5 aswell as D = 6, so the absence of a five-dimensional Nariai instanton on the 
GB branch would seem to suggest otherwise. However, we can't rule out the possibility that 
there are other gravitational instantons whose action is also negative, leading to exponentially 
enhanced production rates. 

To sum up we have demonstrated that the Gauss-Bonnet vacuum is perturbatively un- 
stable in EGB theory. We expect an analogous result to be true generically for all Lovelock 
theories. The result is robust as long as linearised perturbation theory is valid. This is true 
beyond a certain length scale, which grows larger and larger as we approach the Chern- 
Simons limit. Close to that limit, when linearised perturbation theory has broken down 
almost completely, we found strong mixing between the two vacua, so that neither could 
accurately describe the true quantum vacuum. 
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